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1. Introduction 

Given a C°° function, /, on the interval [0, 1] let Rn(/) be the Riemann sum 
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In freshman calculus one learns that 



(1.2) 



What is not perhaps as well known is that if one chooses the U's judiciously, i.e., lets 
fj = j? the O (j?) in l|1.2|l can be replaced by a much better error term, an asymp- 
totic series: 
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in which the Bfc's are the Bernoulli numbers. In particular if / is periodic of 
period 1 the O (4f) in l|1.2|l is actually an 0(N~ ca ). (For an expository account of 
this "Euler-Maclaurin formula for Riemann sums" see |CjS| . ) 

In this article we will prove an n-dimensional version of this result in which the 
interval [0, 1] gets replaced by a convex polytope. We will give a precise formulation 
of our result in fJU however, roughly speaking, it asserts that if A is a simple convex 
polytope whose vertices lie on the lattice, Z™, and if / is in C°°(A) the difference 

J A teJVAnZ" v ' 

can be expanded in an asymptotic series in iV _1 in which the coefficients are ex- 
plicitly computable by recipes resembling . Our formula bears a formal resem- 
blance to the generalized Euler-Maclaurin formulas of |KP| . |KK| . |CSj . |Gu |. |BV| 
et al., however in these so-called "exact" Euler-Maclaurin formulas the functions 
involved are polynomials, not as in the case here, arbitrary C°° functions. Some- 
what closer in spirit to our result is the Euler-Maclaurin formula with remainder 



(1.4) 
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of |KSW| and the Ehrhard theorem for symbols of jGSWj . (Our result also yields 
an Ehrhard theorem for symbols, and its relation to the theorem in |GSW| will be 
discussed in §0|) 

A word about the organization of this paper. In §0we will review the proof of 
the Riemann sum version of Euler-Maclaurin, for the interval, (— oo, 0] and in [21 
show how to extend this result to regions in R n which are defined by systems of k 
linearly- independent inequalities 

(1.5) {ui,x)<Ci, Ui S Z™ , Cj G Z . 

(We will call such regions k — wedges.) 

In [21 we "will derive from this result a Euler-Maclaurin formula for Riemann 
sums over polytopes and in [0 show that our result has an equivalent formulation 
as an Ehrhard theorem for symbols. 

We would like to thank Dan Stroock and Hans Duistermaat for helpful discus- 
sions concerning the material in Section [21 

2. Euler-Maclaurin for the interval (— oo, 0] 
Let t(s) be the Todd function 

In this section we will show that for Schwartz functions, / S S(M.) the difference 

(2.2) ff £,(__) 
has an asymptotic expansion: 

(2.3) ffl + Y(-iy>-i-\f^-V(o)N- 2n . 

V ' 2N ^ y ' 2n)r V ' 

In view of l|2.1|) this formula can be written more succinctly in the form 

("> 5g>(4Hte£) />>*)<*-> 

and it is this version of it which we will prove. 

We first of all observe that if f(x) = e Xx , A > 0, then 
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/ f{x)dx = \e Xh . 



So for N > 2tt\ we may apply the infinite order constant coefficient operator 
- i_ ay 

,N dh/ 



T i~k ~Hh) t° this expression: 



; ' [ 7fWl)Lj ix)dx = T {^Thl a 



T 



N J A 

A e A/ 
N 1 - e- x / N A 



e~£ A e A \ 



\fc=0 



all series being convergent. We conclude that 
(2.5) ±± e -^=(r(±l)f ^(fc-O). 



— OO 



More generally differentiating this identity n times with respect to A we obtain 

M s|B)"«-**= ('(si) (ft = 0) 

verifying i|2.4[l for the function x™e A;l: and hence for the functions of the formp(a;)e Ax 
where p is a polynomial. Now let / be a Schwartz function and p a polynomial 
having the property that f(x) — p(x)e Xx vanishes to order n + 2 at x = 0. Let 



(2.7) g{x) 
Then 



0, x > 

f(x)-p(x)e xk , x<0. 



(2.8) ||3 W (a0lli < oo for i < n + 2 

and by the Poisson summation formula 

( 2 -9) E 4-^)=N 9(Nk). 

— 00<fc<OC — OO<fc<0O 

However, by l|2.8[> 

(2.10) \g(Nk)\ < Const. 7V""fc~ 2 
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for k ^ 0, and 

(2.H) §(0)= f g(x)dx. 

J — oo 

Hence 



(2-12) "TV = / ^)^ + 0(iV-»). 

This shows that l|2.4|l is true for g modulo 0(N~ n ) and hence is true for / 
modulo 0(N-°°). 

Q.E.D 

In we will also need a version of the theorem above for "twisted" Riemann 
sums. Let w ^ 1 be a q th root of unity and let 

1 — uie s 1 — uj 

i>l 

For / S S'(K) we will show that the twisted Riemann sum 

fc=0 

is asymptotic to the series 



(2.14) -J- M + ^^/«( )iV-V 

As above we can rewrite this in the more succinct form 



and we will prove this by essentially the same proof as before: If / = e Xx the 
expression in parentheses is 



(2.16) 



A \ e 



Mi 



1 



N V 1 -, 
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and by setting /i = 0we see that Ij2.15|l is valid for / = e Xx ; and by differentiating 
both sides of (|2.1tj|) by (^r)™ that it's valid for x n e Xx and hence for p(x)e Xx where 
p(x) is a polynomial. Thus, as above, we're reduced to showing that for the function 
g defined by l(2~7jl : 

(2.17) 1 ]T ^(|)=0(iV-"). 

For r = 0, 1, . . . , q — 1, let g T {x) = + Then 



Since 



oo < fc < oo v 7 r— \ — oo<Aj<og 



9r(Nh) = - e % o g 

q V i 

the Poisson summation formula yields, as before, the estimate 



9- 1 poo 

(2.19) Y u \ 9r(x) dx + 0{N~ n ) , 

for the right hand side of l|2.18|l . However, 

/oo 
g (x) dx 
-oo 

and ^'^w 1, — so the first summand in l|2.18|l is zero. Q.E.D. 

We will conclude this discussion of one dimensional Euler-Maclaurin formulas 
by describing analogues of 1|2.4|) and l|2.15|l in which the sum over — oo < k < gets 
replaced by a sum over — oo < k < oo. For simplicity assume that / £ Cq°(1R). We 
claim: 



(2.20) - £ fi^-j = J J(x)dx + 0(N-°°) 

and, for ui a q th root of unity, w^l. 



(2.21) 




= 0(N-°°) . 



To prove (|2. 21)11 we first observe that for c a large positive integer, the left and right 
hand sides of (|2.2UI) are unchanged if one substitutes the function, f(x + c), for /, 
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so without loss of generality we can assume that / is supported on the interval, 
x < 0, in which case (|2.3() is of order 0(N~°°) and (|2.20() is a consequence of l|2.4[) . 
Similarly if we replace f{x) by f(x + cq), with c a large positive integer, the left 
and right hand sides of (|2.21(l are unchanged; so we can assume that / is supported 
on the interval x < 0, and H2.21(l is a consequence of Ij2.15|l . 

3. Euler-Maclaurin for wedges 

Let Z" be the integer lattice in M", (Z n )* its dual lattice in (W 1 )* and (u,x) 
the usual paring of vectors, x € K™, and u £ (R™)*. Given m linearly independent 
vectors, Ui 6 (M™)* we will call the subset of K™ defined by the inequalities 

(3.1) (u^ x) <Ci i = l,...,m 

an integer m-wedge if the Ci's are integers and the u^s primitive lattice vectors in 
{IT 1 )*. Let W be the set (jHU and U the subspace of (W 1 )* spanned by the Uj's. We 
will call W a regular integer m-wedge if ui , . . . , u m is a lattice basis of the lattice 
Un(Z n )* i.e., if 

(3.2) l/n(Z n )*=Bpan z {ui,...,u m }. 
We will need below the following criterion for regularity. 

Lemma 3.1. // holds, Ui, . . . , u m can be extended to a lattice basis, v,\,...,u n 

f(z n y. 

Proof. Let it m +i, . . . , u n be vectors in (Z™)* whose projections onto the quotient of 
(Z n )*by U n (Z)* are a lattice basis of this quotient. □ 

For an integer m-wedge satisfying Ij3.2(l the n-dimensional generalization of Euler- 
Maclaurin is relatively straightforward. 

Theorem 3.2. Let Wh be the subset ofM™ defined by the inequalities 

(3.3) (ui, x) < a + hi , i = l,...,m. 
Then, for / 6 C£°(R n ), 

(«) w S K^)~H^)/ W / W<fa ) (ft = 0) 

where r(si, . . . , s m ) = J]™ i T~(si). 

Proof. By Lemma 13. II we can incorporate u\, . . . , u m in a lattice basis iti, . . . , u n of 
(Z™)*. Let ai,..., a„ be the dual basis of Z™ and let u = YliLi c i a i- Then via the 
map 

(3.5) ieM"^^ 
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one is reduced to proving the theorem for the standard m-wedge: x\ < 0, . . . , x m < 
0, i.e., showing that the sum 

(3-6) ^£/(!>->!) 

summed over all {k\, . . . , k n ) £ Z™, with ki < for i < m, is equal to the expression 



1 d 



h„ 



(3.7) T \Jy~dhJ J "J dxi...dx m J ■J f(x)dx m+1 ---dx n , 

evaluated at ft. = 0, modulo 0(N~°°). Moreover, since the subalgebra of Cq°(R™) 
generated by the products 

f(x) = f 1 (x 1 )...f n (x n ), / 4 eC °°(R), 

is dense in C^°(R") it suffices to prove the theorem for functions of this form, and 
hence it suffices to prove the theorem for n = 1 and m = or 1. However, these 
two cases were dealt with in &\ (See (TO|) and (|!T2U|) .) □ 

We will next describe how 13.4J has to be modified if the condition (|3.2(l isn't 
satisfied. As above let w m +i, ■ • • ,u n be vectors in (Z™)* whose projections onto the 
quotient of (Z ra )* by U (~l(Z™)* are a lattice basis of this quotient lattice. The vectors, 
u\, . . . , u n are now no longer a lattice basis of (Z n )* but they span a sublattice 

(3.8) A* = span z {tti, . ..,it n } 
of (Z™)* of rank n, so the quotient 

(3.9) T = (Z")*/A* 

is a finite group. Let ai, . . . , a n be the basis vectors of M™ dual to ui, . . . , it n . Since 
A* is a sublattice of (Z n )* the dual lattice, 

(3.10) A = span z {a 1 ,...,a n }, 

contains Z™ as a sublattice. Moreover, each element, x £ A, defines a character of 
the group, T, via the pairing 

(3.11) 7 G T -► e 2 ™^ 

and this character is trivial if and only if x is in Z". By a theorem of Frobenius the 
average value of a character of a finite group is zero if the character is non-trivial 
and is one if it is trivial, so we have 



(3.12) Ir £ e ™<T.»> 



1 if x € Z" 
if x 4. Z" . 
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For each 7 € T let 

(3.13) t 7 (si, . . .,s m ) = t Ui (si) ■ ■ -T Wm (s m ) 

where u>k = e 27r ^ 7,afc ^. We will generalize Theorem 13. 21 by showing that for integer 
ra-wedges which don't satisfy condition (|3.2|) one has 

Theorem 3.3. For / G C£° (R n ) 
(3.14) 

^ E / (|) = ^ I) / w /(«) - °> ™ d °(^°°) ■ 

Proof. By (|3.11(l the sum on the left coincides with the sum 

so it suffices to show that the 7-" 1 summand in (|3.14() is equal to the 7- th summand 
in (|3.15|) . Via the map l|3.5[) the 7- th summand in H3.15|l becomes 



(3.16) 




where g(xi, . . . , x n ) = f(v + x\a\ + • • • + x n a n ), and the 7- th summand in (|3.14|) 
becomes 



1 fid 



h± ph 



(3-17) t— ■ r 7 I — — ) / .../ dx x ...dx m / .../ g(x)dx m+ i . . . dx k 



—00 



evaluated at h — 0. (The reason for the factor, l/|r|, is that this is the Jacobian 
determinant of the mapping To prove that l|3.16|l and i|3.17|) are equal 

mod 0(N~°°) it suffices as above to prove this for functions of the form g — 
gi(x) . . . g n (x n ) with gi G Cq°(IR) and hence to show, for i < m 

(3 . 18) ^g„s.(|)~(x„,(iA) £*<».)*.) (».-o) 

and, for i > m 

(3-19) ^E^(ij = ]_Ji{xi)dx i + 0{N-°°), 

and these follow from the identities l|2.15[) and (|2.20() . □ 
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4. RlEMANN SUMS OVER POLYTOPES 

Let A C f be an n-dimensional polytope whose verticies lie on the lattice Z n . 
A is said to be a simple polytope if each codimension k face is the intersection of 
exactly k facets. (It suffices to assume that the vertices of A, i.e., the codimension 
n faces, have this property or, alternatively, that there are exactly n edges of A 
meeting at each vertex.) If the number of facets is d then A can be defined by a 
set of d inequalities 

(4.1) (ui,x) < Ci 

where c; is an integer and Ui € (Z™)* is a primitive lattice vector which is perpen- 
dicular to the i th facet and points " outward" from A. By the simplicity assumption 
each codimension k face of A is the intersection of k facets lying in the hyperplanes 



(4.2) (ui,x) = a ieF 

where F is a k element subset of {1, . . . , d}. Let Wf be the k- wedge 

(4.3) (Ui,x)<a ieF. 

We will say that A is regular if each of these /c-wedges is regular. (As above it 
suffices to assume this for the zero faces, i.e., the vertices of A, or alternatively to 
assume that for every vertex, v, the edges of A which intersect at v lie on rays 

v + tai , < t < oo 

where a\, . . . , a n is a lattice basis of Z n .) 

For regular simple lattice polytopes one has the following Euler-Maclaurin for- 
mula. 

Theorem 4.1. Let Ah be the polytope 

(4.4) (ui,x)<Ci + hi i = l,...,d. 
Then, for f £ C^(R n ) 

(4 - 5) ^A*M<**)L M *>- 0) 

where r(s 1 , ...,s d )= r(si) . . . T(s d ). 

Proof. By a partition of unity argument we can assume that supp/ is contained 
in a small neighborhood of the set (|4.2|l and doesn't intersect the hyperplanes, 
(ui, x) = Ci, i F. Then for i F 

□ 
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However, by (|2.4J) all the terms on the right except the first are integrals of 
derivatives of / over the hyperplane (iii , x) = c% + hi, and hence for hi small arc 
zero. Thus the left hand side of l|4.5[) becomes 



and the theorem above reduces to Theorem 13.21 Q.E.D. 

If A is simple but not regular, one gets a slightly more complicated result. To 
the codimension fc-face of A defined by (|4.2|l attach the subspace 



of (R™)*, the sublattice 



and the finite group 



U F = span R {Mi , i e F} 



Zp = span z {ui , i 6 F} 



T F = U F n (Z")*/Z 



F ■ 



This group coincides with the "torsion group" (|3.9|1 of the wedge W F . Moreover, if 
E is a subset of F, Ue is contained in U F and Z F inZ F , so T F is contained in T F . 
Let be the set of points in T F which are not contained in T F for some proper 
subset, E of F. 

Theorem 4.2. For f 6 Cfi°(M. n ) the sum 

( 4 - 6 ) 4 E / 



N n ^ J \N 

kez n nNA 



k 



is equal mod 0(N °°) to 



i>; >>7i^) / 1 1'' = "»- 




Proof. As above it suffices to prove this for supp/ contained in a small neighborhood 
of the set (|4.2p . and not intersecting the hyperplanes, (ui,x) = Ci, i ^ F. Then as 
above, the only contribution to the sum (|4.7|l is 



j>(^)L/H , ''= 0) 



and Theorem 14 . 21 reduces to Theorem 13.31 □ 
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5. An Ehrhart theorem for symbols 

A function, / £ C°°(R™) is a polyhomogeneous symbol of degree d if, for large 
values of x, it admits an asymptotic expansion 

— oo 

(5-1) 

whose summands are homogeneous functions fj € C°°(M™ — {0}) of degree j. Let 
/ be such a function and let A be a simple lattice polytope in K n containing the 
origin in its interior. The Ehrhart function of the pair, /, A, is defined to be the 
function 

E(f,A,N)= ]T f(k), N E I*+ . 
1-eAfAnZ" 

In |TTSW it was shown that 

E(f,A,N)- f fdx 
Jna 

had an asymptotic expansion 

— oo 

(5.2) c J Nj + c 

j—n+d 

for N large. 

The main result of this section is a variant of this result. As above let A be 
a simple lattice polytope in W 1 and let Ca be the polyhedral cone consisting of 
all points, (x±, . . . , x n , x n+ i), in M n+1 with x n +i > and (xx, . . . , x n )/x n +x G A. 
Then, for N £ Z+, A^A is just the slice of Ca by the hyperplane, x n +x = N. We 
will prove that if / € C°°(R" +1 ) is a homogeneous symbol of degree d the sum 

(5.3) J2 /( fc ) 
admits an asymptotic expression of the form H5.2|l . 

Remarks. (1) This result, albeit very close in spirit to the theorem in _GSW 
cited above, doesn't, as far as we can see, seem to be a trivial consequence 
of it. 

(2) This result has a number of applications to spectral theory on toric varieties 
which we'll explore in future publications. 

(3) As a corollary of this result one gets another variant of the Ehrhart theorem: 
Let 

A" := {(xx,---,x n+1 ) e C A , x n+1 < 1}. 
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This (n + l)-dimcnsional polytope is not in general simple. However a 
version of theorem described at the beginning of this section is still true, 
namely 

— oo 
i— rf+n+1 

as one can see by summing the differences 

E(f,A*,N)-E(f, A»,JV-1) 

and noting that each difference is exactly Ij5.3|) . By combining this result 
with the Danilov "desingularization trick" |Daj one can extend the Ehrhart 
theorem to a much larger class of convex lattice poly topes. We will discuss 
the details elsewhere. 

Proof. As above let 

— oo 
i—d 

where fi(x\ 7 . . . , x n +i) is a homogeneous function of degree i. Then on the cone, 

r / \ iff ^1 i 

/H^l) • ■ ■ j #n+l ) — x n-\-l Ji I ) * * " i ) * 



so if we set f%{x\^ . . . , x n ) — /i(xi, . . . , x n , 1) the sum (|5.3|) is equal to the sum 

( 5 - 4 ) Ni E 

which is N' l+n times the Riemann sum 

Thus, by Theorem l4.2l each of these summands admits an asymptotic expansion: 



E c *.^ fc 

and hence so does the sum (|5.3|) . □ 
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